Spontaneously Symmetry-Broken Current in Coupled Nanomechanical Shuttles 
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We investigate the transport and the dynamical properties of tunnel-coupled double charge shut- 
tles. The oscillation frequencies of two shuttles are mode-locked to integer multiples of the applied 
voltage frequency uj. We show that left/right-symmetric double shuttles may generate direct net 
current due to bistable motions caused by parametric instability. The symmetry-broken direct cur- 
rent appears near uj — flo/(2j — 1), (j = 1,2,.. .), where fio is the dressed resonance frequency of 
the relative motion of the two shuttles. 
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Recently, nanoelectromechanical systems (NEMS) 
have attracted great attention due to fundamental as- 
pects in new electrical transport properties as well as new 
technology complementary to conventional MEMS engi- 
neering jlj. A prototype of NEMS, single electron shut- 
tle, was suggested by Gorelik et al. Q which is a single 
electron transistor combined with its mechanical degree 
of freedom. The charge shuttles can be realized in two 
different ways in experiments. In a top-down approach, 
the charge shuttles can be realized by miniaturization 
of silicon structure 0, 0] . The other way is a bottom- 
up approach, where the charge shuttle is produced from 
molecules such as Cqq 0- 

In this Letter, we report a theoretical study on the 
transport properties and the dynamics when two charge 
shuttles are coupled through electron tunneling. The 
nonlinearity involved in this system shows an interest- 
ing bistable regime, where the electric current character- 
istics are of potentially great importance in NEMS ap- 
plications. An interesting situation arises when the two 
shuttles are totally symmetric. In this case, through a 
dynamical symmetry breaking the system produces a net 
direct electric current. Symmetry- broken electric current 
under time-periodic perturbation has been an interesting 
topic of many theories and experiments 6] where a com- 
plete understanding is still necessary. We show that the 
double charge shuttle (DCS) allows for symmetry-broken 
DC current caused by parametric instability. This is in 
contrast to the single shuttle system where the DC cur- 
rent is vanishing in symmetric configuration. 

We start with the formalism of single electron tunnel- 
ing and classical dynamics which has been used for sin- 
gle charge shuttle The drain lead is in the left side 
of the source lead where the ac voltage V(t) = Vq sin (u)t) 
is applied to the source compared to the drain, see Fig.^ 
The capacitance is supposed to be not sensitive to the 
displacement while the resistance is a function of dis- 
placements x\ and x 2 : 
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FIG. 1: A schematic hgure of the double charge shuttle. 
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Here, A is a phenomenological tunneling length. When 
the mutual capacitance of the j-th junction is Cj, 
the capacitance matrix is constructed as Cm = + 
Ck+i, — Cfe+i, — Cfc, for I = k, I = k + 1, 1 = k — 1, other- 
wise, respectively. The internal charging energy is given 
by e(Q 1 ,Q 2 ) = \ Y?k,i(. c ~ 1 )kiQkQi- The energy loss Ej 
of the j-th junction is E 1 = %&eV + e(Qi, Q 2 ) - e(Qi - 
), E 2 = ^eV + e{Q x , Q 2 ) - e{Q x + e,Q 2 ~ e), and 
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E 3 = J£±eV + e(Qx,Q 2 ) - e{Qi,Q 2 + e) with the total 
capacitance ctot- The rate of the tunneling from left to 
the right at the j-th junction at low temperature and 
environmental impedance is 
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where the island charges Qi,Q 2 are integer multiple 
of — e; Qi = —rue. By replacing V,Qi,Q 2 with 
—V, —Qi, —Q 2 in the right-hand side of the above equa- 
tion, one can get *Tn x ,n 2 - 

Now we come to the equations of motion. The 
time evolution of the probability P ni ,n 2 f° r the island 
charges Q\ = —n\e, Q 2 = — n 2 e is given by the follow- 
ing rate equation 
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FIG. 2: (a) The symmetry-broken DC current in units of 7o = 
Vb /R vs. frequency of the applied AC voltage. The solid line 
refers to the result from the full calculations mentioned in the 
text and the dashed line refers to the result from the simplified 
calculation in the adiabatic limit. The parameters are fixed to 
c = 100 aF, 7? = h/e 2 , A = 20 nm, L = 500 nm, lu = 50 MHz, 
7 = 0.025 w , m = 1.1 x 10" 21 Kg, and V = 10.8 mV. (b) The 
corresponding stroboscopic plot of x — xi ~ 22 obtained for 
a fixed phase of the oscillating voltage. 
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At low temperatures the fluctuation of the displace- 
ments X\,X2 are negligible compared to the charge fluc- 
tuations. While the charge fluctuations are important to 
the noise properties, we suppose 1/T be much smaller 
than the typical time scale of displacement and inves- 
tigate the mechanical motions using the mean island 
charge (Q,-(i)) = -e(n,(t)) = -eJ2 ni ,n 2 niP(ni,n 2 ,t). 
The island charges experience the force produced by the 
electric field -V(t)/L; 
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where L is the source-drain distance and 7,, ujQi denote 
the friction constant and the natural angular frequency 
of i-th shuttle. Note that the resonators are modelled 
as linear oscillators. The nonlinearity of the full system 
comes from the coupling via tunneling. The electric cur- 
rent from the source in the right-hand side to the drain 
in the left-hand side is computed as 
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In the following we will consider a symmetric config- 
uration: LU 0i = LU 0l 7i = 7, TOj = TO, i?l(0) = i? 3 (0) = 

0.5i? 2 (0) = R, and c x = c 3 = 2c 2 = c. In Fig. 12a), we 

plot the time-averaged current Idc — 57 Iq^^ I(i)dt. 
From the symmetry of the configuration one would ex- 
pect that the DC current computed according to Eq. © 
is zero for all frequencies of the applied AC voltage. 
Fig.|2ta), however, shows that in the vicinity of the natu- 
ral frequency of the oscillators, u>q, the DC current can be 



negative or positive depending on the initial conditions. 
The initial conditions are specified by the initial position 
xi, X2 and the velocities x\, £2 at a given starting time. 
In the regime of zero current all initial conditions are 
attracted to an unique stable asymptotic solution (the 
so-called attractor). In the nonzero current case, two 
such attractors coexist and cause a dynamical symmetry 
breaking. 

The global behaviour of the attractors can be conve- 
niently represented on a stroboscopic section (a special 
case of a Poincare section 0) defined by a fixed value of 
the phase of the driving voltage. Fig.|2{b) shows the stro- 
boscopic section at a fixed phase of the applied voltage. 
We find the center of mass coordinate does not move; 
X = (xi + x-i)j1 » for all frequencies u>. Meanwhile in 
the regime of nonzero DC current we find two symmetry- 
related solutions for the relative coordinate x — x± — a?2. 

We may get further insight by investigating the adia- 
batic limit of the system. In Fig.[3we plot the DC current 
and relative coordinates computed from the adiabatic ap- 
proach using dotted lines. One can notice that the full 
numerical results are quite close to those obtained in the 
adiabatic limit. Therefore we may rely on the adiabatic 
approach for the analysis of the bistability. 

In the adiabatic limit the electronic relaxation is much 
faster than the mechanical motion, i.e. 

ojRc , ujqRc < 1. (5) 

In this case, a classical circuit analysis of Fig. 2] on the 
charges qj accumulated in each capacitor Cj gives 
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From the solutions of the above linear equations, the net 
charges of each islands Q\ = 91—92, Q2 = 92 — 93 are 
simply given by 
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Plugging these into the equation of motion for the me- 
chanical degree of freedom © gives the equation of mo- 
tion for the center of mass coordinate 
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where x = x\—X2 is the relative coordinate. Clearly, X = 
is a solution, in agreement with our numerical finding 
that the center of mass is not moving. Exploiting this 
fact, the equation of motion for the relative coordinate x 
can be derived as 



x + jx + oj q x 
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From the above equation one can see that x = is a 
trivial solution. However, this solution can be unstable 
as we will prove in the following. Note that if a nontrivial 
solution of x{t) exists then — x{t) is also a solution. This 
is the pair of bistable solutions. Since the equation is in- 
variant under the time-translation operation t —> t+n/uj, 
a periodic solution should satisfy x(t + ir/u>) — ±x(t). 
As will be shown later, this parity is important for the 
nonzero DC current. 

The electric current in the adiabatic limit can be 
expressed as I(t) = qi(t)/[i?i(x)ci] = V(t)/[R x (x) + 
R2{xi,x 2 ) + ^3(^2)]- Therefore the time-averaged DC 
current reads 
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where Iq — Vq/R. The dashed line in Fig. EJa) was ob- 
tained using the above formula which also shows a clear 
bistability. 

Now we are going to show that the origin of the 
symmetry-broken current and the bistability is paramet- 
ric instability. By linearizing the term tanh(3x/4A) in 
the right hand side of Eq. 1)100. we get 
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frequency of the relative motion; 



and £!n is the dressed harmonic 
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Equation (|12f> is called a damped Mathieu equation 
which is a paradigm for studying parametric reso- 
nance The stability analysis in Ref. 0| up to sec- 
ond order of M shows that the motion is unstable in the 
interval c<j_ < u> < lu + where 
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with M 



This interval is called principal in- 
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stability interval. Obviously, the interval has finite width 
only if M > 7/^0, he. the strength of the driving, Vo, 
must be sufficiently large in order to get parametric in- 
stability. 

From the above analysis and Fig.EIb) we conclude that 
when entering the principal instability interval from the 
right hand side the unique attractor turns into an unsta- 
ble solution (a repeller) thereby creating two new attrac- 
tors. This scenario leading to bistability is termed a su- 
percritical pitchfork bifurcation. The bifurcation occur- 
ing when entering the principal instability interval from 
the left hand side is different which can be seen from 
the abrupt change of the value of x. Here we have three 
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FIG. 3: Bifurcation diagram /j, vs. uj/ujq in the adiabatic 
limit for 7 = 0.025ljo- Bistable regions (Arnol'd tongues) 
are mode-locked with p = 1, 2, 3, 4, . . .. The nonzero DC cur- 
rent arises only when p is odd. The dotted line denotes the 
boundary of the principal instability region obtained from the 
analytic formula 1141 . 



attractors and two repellers (not shown). The repellers 
and the central attractor [x = 0) merge in a subcritical 
pitchfork bifurcation creating a repeller at x = 0. 

In Fig. El we plot the phase diagram of /i and lv/ujo 
which indicates the locations of the bistable regimes com- 
puted from Eq. I|1(J|I . The dotted line denote the an- 
alytic results in Eq. i|14[l from the linearized equation 
(|12fl . The reasonable agreement of the parameters show- 
ing the instability confirms that the parametric amplifica- 
tion causes the instability in the system. Beside the prin- 
cipal instability region there are higher-order regions in 
which the ratio of the frequency of the oscillators and the 
applied voltage is mode-locked to p = 2, 3, These re- 
gions, and also the principal instability region, are called 
Arnol'd tongues |ll| . 

The second interval of instability denoted by p = 2 in 
Fig. |3 arise from second harmonics of the system. The 
analytic formula for this region from the stationary solu- 
tions in Ref. 01 g ives 
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The higher order bistable regions are found in our nu- 
merical calculations. The center of the intervals in each 
bistable regions are located where 
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The main characteristics of the each bistable region is 
that the (unharmonic) frequency of the relative motion 
is given by pu> even though there is a mismatch between 
the natural harmonic frequency and the frequency of the 
voltage, which is called mode-locking in the field of non- 
linear dynamics 0. In Fig. H we plot the time series 
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FIG. 4: Time evolution of the relative coordinate x = x\ — X2 
for typical sets of values (w/wo> /-0 belonging to the bistable 
regimes of p = 1,2,3,4; uj/uj = 1.45,0.79,0.53,0.415 and 
fi = 1.414. 



of the relative coordinate x for various frequencies corre- 
sponding to different bistable regimes. One can clearly 
see the p-th order mode-locking. For instance, in spite 
of the large mismatch between u> and loq when p = 1 
(to = 1.45wo) the oscillation frequency is clearly given 
by u. 

It is important to note that the rectified DC current 
arises only when p is an odd number of integer. When 
p is even, the periodic solution x(t) has even parity i.e. 
x(t+Tr/uj) = x(t). In this case the integration in Eq. Hll|) 
from to to t + tt/lo is cancelled by the integration from 
<o + 7t/cj to to + 2tt/uj. 

One may ask whether the parametric instability can 
cause symmetry breaking already in a single shuttle sys- 
tem. When we write the tunnel resistances as Ri(x) — 
R(0)e +x / x and R 2 (x) = i?(0)e~ x / A where x is the coor- 
dinate of the shuttle, the equation of motion is given by 

x + jx + luqX = ~ cV °mL ~ tanh f similar to Eq. (|10|l . 
Therefore, in the sense of the mechanical motion, the 
parametric resonance gives rise to the bistability and the 
bifurcation diagram is quite similar to Fig. However, 
the single shuttle systems does not allow a symmetry- 
broken DC current. The absence of the DC current in the 
symmetric single shuttle is clear from the time-averaged 
current formula of the system 
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Note that the symmetry of the periodic solution x(t + 
tt/lu) — ±x(t) always ensures cosh [x(t + n/ui)/X] = 
cosh [x(i)/A] confirming that the above integral is zero. 

In experimental systems, the charge and displacement 
fluctuations exist due to the discreetness of charges and 
finiteness of temperature. Noise properties of two col- 
loidal particles in DC source-drain bias have been nu- 
merically investigated in Ref.01 where the shot noise 
produces random telegraph noise. The noise properties 



in the presence of oscillating source-drain voltage are ex- 
pected to contain interesting information on the mechan- 
ical properties, which is currently under study. For in- 
stance, the shuttles are expected to show enhanced noise 
power in the bistable regimes we studied in this work. 

Our results are robust against the finite temperature 
effect which modifies the tunneling rate in Eq.^.This 
is because the dynamical symmetry breaking appears as 
a classical effect. While the perfect left/right symme- 
try does not exist when the tunnel resistances differ by 
multiple of 5 or 10, we find the effects are still visible 
from the enhancement of the net current. The enhanced 
net current is also visible when the resonance frequen- 
cies and tunneling coefficients differ by 10 percents. This 
robustness comes from the nonlinear phenomenon, mode- 
locking, where the mechanical motions are locked to the 
driving voltage in spite of small variations of the system 
parameters. One may also note that the tunneling length 
A enters the instability condition through (i, so the net 
current can exist in a reasonable range of A if /J, and u 
belongs to the Arnol'd tongue. 

In summary, the transport through two tunnel- coupled 
symmetric charge shuttles has been studied. We found 
that the oscillation frequencies of the two shuttles are 
mode-locked to the applied oscillating voltage. Moreover, 
we observed a dynamical bistability which allows for non- 
zero electric currents. The origin of this phenomenon has 
been traced back to parametric instability induced by the 
nonlinear coupling of the mechanical and electrical degree 
of freedom. 
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